Collective Cooper-pair transport in the insulating state of Josephson junction arrays 
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We investigate collective Cooper-pair transport of one- and two-dimensional Josephson junction 
arrays in the insulating state. We derive an analytical expression for the current-voltage character- 
istic revealing thermally activated conductivity at small voltages and threshold voltage depinning. 
The activation energy and the related depinning voltage represent a dynamic Coulomb barrier for 
collective charge transfer over the whole system and scale with the system size. We show that 
both quantities are non-monotonic functions of magnetic field. We propose that formation of the 
dynamic Coulomb barrier as well as the size scaling of the activation energy and the depinning 
threshold voltage, are consequences of the mutual phase synchronization. We apply the results for 
interpretation of experimental data in disordered films near the superconductor-insulator transition. 
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Recent experimental studies of the superconductor- 
insulator transition (SIT) (see (ij for a review) in thin 
disordered superconducting films proved formation of a 
collective insulating state exhibiting thermally activated 
Arrhenius-like conductivity at low biases [1, |3, 0] and 
the threshold voltage depinning [3, @| behavior. Be- 
low the depinning voltage, Vt, a film falls into a zero- 
conductivity phase and abruptly switches to a finite con- 
ductance regime when the bias achieves Vt ■ The discov- 
ery of the novel phase existing in a narrow window of 
disorder strength near SIT poses a quest for a compre- 
hensive theory. 

Josephson junction arrays serve as a perfect testing 
; ground for SIT studies (see, e.g., 0, 0, E H H)- A 
salient similarity of the voltage threshold behavior in SC 
films [1, 0] and the voltage depinning in one-dimensional 
Josephson arrays suggests an intimate relation be- 
tween these systems. Further parallel appears from the 
striking observations of voltage threshold Vt dependence 
on the array length in [lot . the saniple size dependent 
activation energy, ksTo, observed in M, and the connec- 
tion between Vt and Tq revealed in [3|. An advantage 
of a Josephson array as a model system is that it offers 
a straightforward theoretical description of the current- 
voltage characteristics, which is what measured in all the 
major experimental studies of SIT. In this Letter we de- 
, velop a theory of the collective transport of large Joseph- 
son junction arrays in the insulating state and apply our 
results for interpretation of experimental data on SIT. 

The current-voltage characteristics of Josephson sys- 
tems in an insulating state were discussed in a single 
junction ll,[l3,[I3| and two junction 14, [l^ systems. 



Each junction is characterized by the Josephson cou- 
pling energy, Ejq = hlc/2e, where Ic is the Josephson 
critical current, and by charging energies Ec related to 
inter- island capacitance and Eco associated with capac- 
itance to ground, Co [16|. We consider an insulating 
state, where charging energies Ec,Eco » Ej. We as- 
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FIG. 1: (Color online) Sketch of the considered array ge- 
ometries. An external current I is injected from the left 
through the electrode having the superconducting phase X-c 
and extracted through the right electrode with the phase 
Xr- Upper panel: One-dimensional array of superconduct- 
ing islands (squares) connected by two Josephson junctions 
(crosses) to neighbors corresponding to experimental system 
of Lower panel: Two-dimensional M x N Josephson 

junction array. 



sume further the superconducting gap A > Ec- This 
implies that the transport is mediated by the thermally 
activated motion of the Cooper pairs 



17| . We show that 



in the regular and disordered arrays a collective current 
state develops. This state is characterized by the energy 
gap, Ac, stemming from the coherent Coulomb blockade 
effect, i.e. the Coulomb blockade involving all junctions 
and extending over the whole system. We derive a low 
bias I-V dependence: 
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Eq.iJT]) reveals that there are two dynamic regimes: first. 
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thermally activated charge transfer with the resistance 



R oc exp[AJ{2kBT)] 



(2) 



at eV ^ Ac and, second, threshold behavior at V ^ 
Vt — Ac, where the activated conductivity turns to 
a finite non-activated transport. We find that in a 
regular ID array Ac ~ EcL/d, while in a 2D array 
Ac ~ Ec \n{L/d), where L is the array length and d is the 
size of the elemental cell of the array. We show that at 
the depinning voltage threshold in disordered 2D systems 
the current breaks through the system along the first per- 
colative path connecting the leads. As a result, the de- 
pinning transition acquires an one-dimensional character, 
and the depinning voltage Vr, scaling linearly with the 
sample length, much exceeds the corresponding low bias 
activation energy. We demonstrate, finally, that the mag- 
netic field perpendicular to the film gives rise to a non- 
monotonic Ac{B) dependence in an excellent agreement 
with the experimental data on both ID artificial Joseph- 
son arrays and superconducting films near SIT [1, . 

Let us consider N x M superconducting islands com- 
prising a two-dimensional array closed by a small (as 
compared to the quantum resistance for Cooper pairs 
Rcp = h/Ae^ ~ 6.45 kfi) external resistance, Rext^ see 
Fig. I. We assign the fluctuating order parameter phase 
Xijit) to the {i,j}-th superconducting island (see Fig. 
I). The phases of the left- and right leads, Xi(^) and 
Xnit), respectively, are fixed by the dc voltage V across 
the array: 



XR-XL^2eVt/h + ^{t) , 



(3) 



where ip{t) describes fluctuations in the leads. We single 
out the leftmost, « = I, and rightmost, i = N , columns of 
islands directly coupled to leads and represent the array 
Hamiltonian in a form: 
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Here 



^0= E 



r 



iXv - Xkif ~ Ej cos(xij - Xki) 



the brackets (ij, kl) denote summation over the pairs of 
adjacent junctions, and the last term in ([5]) represents 
the self-charge energies of superconducting islands. The 



Hint term in (4) describes coupling of phases on the leads 
to the thermal heat bath 11 1 . 

The dc Josephson current through the array is 



1 r 

Is{V) = Ic lim - dty 



Xij (t) +XNj{t) 



X sm 
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where the brackets (...) stand for an averaging over ther- 
mal fluctuations in the leads and quantum mechanical av- 
eraging over phases of internal junctions {xij (t)) and over 
the variable 4>j = {xij +XNj)/'2- We construct the time- 
dependent perturbation theory with respect to small pa- 
rameter Ej/Ec, similarly to the case of the Cooper pair 
two-junctions transistor [l3,[i3li omitting the last term 
in ([5]) since in most experiments C ^ Co, and thus 
Ec ^ Eco- In the first order one finds: 
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2eVt/h + yj{t) + xi{t) - XN{t) 



In the second order, using the approach developed 
one arrives at: 
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IsiV)=MI, 
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where 5 = Ae^RextkBT refiects the Gaussian character of 
the current noise in the leads due to thermal fluctuations 



111 [18|. The correlation function of internal phases is 



defined as 

K{t) = {exp{i[xij{t)-XijiO)~XNj{t)+XNjm))Ho- (9) 

In the two-junction system (single Cooper pair transis- 
tor), xii = XNji K{t) = 1, and we recover the re- 
sults of 15]. In the zero-approximation one neglects the 
Josephson coupling inside the array, and K{t) can be 
found in a closed form as an analytical continuation of 
K{t), where r is the imaginary time: 
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Expanding phases Xiji^) over the Matsubara frequen- 
cies Urn = 2'KkBTm/h as Xiji^) = X] exp(ia;„f)xy (t^m)) 
and going over to charge representation, Xiji^m) — 
xij 1 (5) nij(2£;cfc_Br/(fttj„i)2)[exp(-iLj„iT) - 1], with being 
the number of Cooper pairs localized on the {zj}-th 
island, one eventually obtains the correlation function 
K(t): 



K{t) = exp ( - 2AckBTt'^/h'^ - 2iAct/h) 



(11) 



(10) 
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where Ac is the barrier for the Cooper pair propagation 
through the whole system defined by the relation 
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Importantly, when deriving Eqs. (|10j|12p . the non-zero 
windings numbers Wij = [xijih/ksT) - Xy(0)]/(27r) 
were neglected. This holds only at not very low 
temperatures, T > Ec/iirkB)- At T = Ec/{iTkB) 
a Berezinskii-Kosterlitz-Thouless-like transition into 
a super-insulating state with the resistivity R oc 
exY>{{Ac/ Ec) exp[Ec/ (nkBT)]} occurs In what fol- 
lows we consider the interval of moderately low temper- 
atures, EcKukB) <T < Ac/(2fcs). 

Plugging pT|) into ([8]), one finally derives the current- 
voltage characteristics of the insulating state as given by 
Eq.([T]). At low biases, V <^ Ac, it gives an Arrhenius ac- 
tivation temperature dependence ^ for the resistance, 
with the activation energy, Ac, defined by Eq. Car- 
rying out calculations one finds for a regular array: 

^ _ UEc/2)min{Xc,L}/d, for ID arrays, ^^^^ 
^ I (_Ec/7r) ln(min {Ac, for 2D arrays. 



where Ac — dy^Eco/Ec is the screening length related to 
capacitance to ground [20| . Usually, in experiments Eco 
is so high that Ac well exceeds the sample size L. 

Dependence of the activation energy on the sample size 
was recently observed in InO superconducting films 
Plotted in Fig. [2^ are activation energies, Tg = Ac/fcs, 
extracted from Fig. 4a of [5, vs. logL. One clearly sees 
logarithmic scaling of Tq as Eq. (fTS]) predicts. 

The nature of the Coulomb barrier Ac and its size scal- 
ing can be understood in terms of the mutual phase lock- 
ing or phase synchronization in the Josephson junction 
array. In the Coulomb blockade regime the charge at each 
junction is fixed, and, therefore conjugated phases fluctu- 
ate freely. Yet, the exponentially small dc Josephson cur- 
rent couples phases of the adjacent junctions to provide a 
minimal power dissipation in the array. This establishes 
a global phase-coherent state, and transport occurs as 
a co-tunneling of Cooper pairs through the whole array. 
The probability of such a process in an ID array is pro- 
portional to [exp(— E'c/fcsT)]^, giving the total Coulomb 
barrier as Ac ~ EcN. Another way of thinking is to 
say that synchronization builds on the large screening 
length Ac which allows for the small charge fluctuations 
at each junction to interact over the whole system. In 
this sense the linear (logarithmic) scaling of Ac reflects 
linear (logarithmic) growth of the Coulomb energy in ID 
(2D, respectively) systems. As a result, synchronization 
is rigid with respect to disorder: even large (of the order 



of the quantity itself, but Gaussian) fluctuations in Ec, 
Eco, and Ej, as well as the effect of the offset charges 
is negligible as compared to the huge magnitude of Ac. 
That is why this scaling of Ac holds even in the amor- 
phous superconducting fllms |4|, where the granularity 
is of a self- induced nature 0,123] and the variations in 
Josephson coupling strength are small. 

The current-voltage characteristic of Eq. ([T]) is valid as 
long as (Ac — eV)^ ^ 2AckBT. At temperatures of in- 
terest, T < Ac/ {2kB), this gives an accurate estimate for 
the threshold voltage of the regular Josephson junction 
array as: 

eVr 0^ Ac , (14) 

with Ac from ([T3|) . This result holds in disordered sys- 
tems as well. However, in both ID and 2D systems the 
threshold voltage scales linearly with the size of the sam- 
ple, i.e. the ID scenario works. Consider first an ID 
chain. The current state forms at the threshold voltage 
as a result of the dielectric breakdown where the collec- 
tive charge transfer over the whole array occurs. The 
associated Coulomb energy Ac ^ {L/d){Ec) again scales 
linearly with the system size as long as the distribution 
in Ec is not exponentially broad and the average {Ec) is 
well defined. The size dependence of Vt on the sample 
size was observed in [13], where the chain of SQUIDs, 
schematically shown in the upper panel of Fig. 1, was 
studied. One can see from the Fig. 2d of |lfl] that for 
two largest samples indeed Vr oc L. In 2D arrays with 
disorder the dielectric breakdown becomes of a percola- 
tive nature and occurs along the first "lowest resistance" 
path connecting the leads. This retains a ID scaling of 
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FIG. 2: (Color online) (a) Activation energy To = Ac/fcs 
plotted as function of logarithm of the sample length (squares 
are the data from 4j). (b) The experimental data from Q: 
activation energy (circles, left axis) and voltage threshold 
(squares, right axis) as functions of the magnetic field B. 
The lines are are according to Eq. (|15|l : Ej^{B), with values 
aEjo/Ec = 0.8 and Aioop = 1.4-10-^ fim'^ fits To (solid line); 
Vt{B) is the ID quantity and is fitted by -Ej^ with the same 
Aioop and aEjo/Ec = 0.96 (dashed line) reflecting slightly 
different geometric factor. 
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Vt- Consequently, one expects that in 2D films the corre- 
sponding energy eVr is much larger than the activation 
energy determined from the low bias resistance behav- 
ior (12) ■ Indeed observed in [sl was eVr/AsTo w 220 at 
the magnetic field 0.7 T. The above percolative picture is 
identical to the threshold charge depinning in 2D arrays 
of metallic dots investigated numerically in ^2j\. 

Next, we discuss the effect of the magnetic field on the 
activation energy and voltage depinning threshold. The 
field modulates the effective Josephson coupling: in ID 
SQUID chain one has Ej^ = Ejo \ cos(7r/)|, while in the 
2D array i;}^ = Ejo{l-Afsin^[7r{l-f)/A]} [H, where 
/ = eBAioop/i^h, Aioop is the area of either the elemental 
SQUID or the plaquette in the 2D array. The correction 
to Coulomb barrier in the first order perturbation theory 
with respect to Ej/Ec follows from (|llll2p : 



^,{B)^^,[l-aE.,{B)/E, 



(15) 



where the parameter a is of the order of unity, and de- 
pends on the geometry of the lattice. The field modula- 
tion of Ej{B) yield non-monotonic field behaviors of Tq 
and Vt{B). Shown in Fig. 2b are fits to activation energy 
Tq and Vt vs B dependencies to the experimental data 
from 0. The quantity aEjo/Ec = 0.8 is chosen to match 
To(0) to i? = experimental value and reflects that ex- 
periments were carried out in the vicinity of SIT (still 
allowing "borderline" estimates within the perturbation 
theory). The loop area is defined unambiguously by the 
position of the maximum in Tq (only the branch corre- 



sponding < / < 1/2 should be taken in E"^^ [23). Us- 
ing the same Aioop, the theoretical VriB) matches with 
the data of Q (Fig. 2b); this procedure means the effec- 
tive absence of fitting parameters. Fig. 2 confirms the 
2D nature of activation energy and the ID scenario of 
depinning threshold. 

In conclusion, we have developed a theory of collec- 
tive Cooper pair transport in the insulating state of one- 
and two-dimensional Josephson junction arrays. We have 
obtained the Arrhenius low-bias resistance and derived 
the corresponding activation energy. We have shown 
that both, the activation energy and the voltage depin- 
ning threshold, represent the dynamic Coulomb barrier 
Ac controlling collective charge transfer in the insulat- 
ing state. In Josephson junction chains the activation 
energy and voltage threshold coincide and both scale lin- 
early with the chain length. In two-dimensional arrays 
the activation energy scales logarithmically with the sam- 
ple length, while threshold voltage, Vt, exhibits the ID 
linear scaling, since disorder sets the percolative dielec- 
tric breakdown mechanism of charge depinning. We have 
proposed that the physical origin of the energy gap and 
its scaling is the mutual phase-locking in junction arrays 
which maintains even in disordered systems. We expect 
that at temperatures above the energy gap Ac/fcs the co- 
herent synchronized state breaks down and the collective 
activation transport transforms into a usual local variable 



range hopping as observed in [1, [S^l . We have demon- 
strated that modulating Josephson coupling by the mag- 
netic field leads to a peaked VT(i3)-dependence in agree- 
ment with the experimental findings for ID Josephson 
arrays (lo| and for superconducting films near 
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